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^ ■ We propose a new nucleosynthesis process in which various nuclei are formed 

through the liquid-gas phase transition of supernova matter as a preprocess of the 
standard r-process, and study its possibility of realization qualitatively. In the 
relativistic mean field and statistical model calculations, we find that this process 
may take place and help the later r-process to proceed. 
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\Q ' §1. Introduction 

(N ■ 

i It is generally believed that there exist several phases of nuclear matter. Among 

' the phase transitions between them, the nuclear liquid-gas phase transition has been 

CN ■ extensively studied in these three decades. It takes place at relatively cold (Tc = 

(5—8) MeV) and thin [ps ~ Po/S) nuclear matten and it causes multifragmetation of 
expanding nuclear matter in heavy-ion collisions liP . In the universe, the temperature 
and density region of this transition would be probed during supernova explosion. 
O I Since supernova explosion can be regarded as expansion of infinite matter composed 

of nucleons and leptons, its fragmentation through the phase transition may produce 
various nuclei including heavy nuclei which are believed to be synthesized in r-process 
or rp-process. Then the phase transition during supernova explosion may be helpful 
^ ' to solve some problems in these standard processes. 

^ . In this work, we propose a new kind of nuclear synthesis process through the 

liquid-gas phase transition (LG process) as a preprocess of the usual r-process. In the 
early stage of explosion, since the density and the temperature are high enough to 
keep statistical equilibrium and to trap neutrinos, entropy per baryon S/B and lepton 
to baryon ratio Yl are to be conserved. This condition is expected to be satisfied 
within the neutrino sphere {pB = (10~^ — l)po)- At these densities, supernova matter 
expands almost adiabatically and cools down. If the critical temperature Tc of the 
liquid-gas phase transition is higher than the freeze-out temperature, this matter will 
fragment and form various nuclei in a critical manner through the phase transition. 
As the baryon density decreases, charged particle reactions are hindered then the 
chemical equilibirum ceases to be kept, namely the system freezes out at this point. 
The statistical distribution of fragments at freeze out will give the initial condition 
of the later processes such as the r-process. 
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Although the importance of the liquid-gas phase transition in supernova mat- 
ter was already noticed and extensively studied before the main interest in these 
works was limited to the modification of the equation of state (EoS), then the nuclear 
distribution as an initial condition of the r-process was not studied extensively. In 
addition, the mean field treatment assuming one dominant fragment configuration 
was applied rather than to consider the statistical ensemble of various fragment con- 
figurations, which is important at around the critical temperature for sub-saturation 
densities. 

In this study, we first investigate the liquid-gas coexistence region in the {pB,T) 
diagram and the ejection possibility of materials inside this coexistence region within 
the Relativistic Mean Field (RMF) theory. Next, by using a statistical model, we 
evaluate the fragment formation at around the phase transition and in the coexistence 
region. Then we compare the fragment distributions at various {pB,T) points with 
the solar abundance. 

§2. Relativistic Mean Field 

Here, we study the phase diagram of 
supernova matter by using RMF theory. 
For simplicity, we have made following 
three assumptions: (1) Supernova mat- 
ter is composed of nucleons, electrons, 
electron-neutrinos and their antiparti- 
cles. (2) Neutrinos are trapped then the 
lepton to baryon ratio is conserved. (3) 
Two coexisting (liquid, gas) phases are 
of infinite size. Based on these assump- 
tions, we have carried out the RMF cal- 
culation with the interaction parameter 
set TMIQ), which reproduces the bind- 
ing energies of stable as well as unsta- 
ble nuclei. In supernova matter, there 
are three conserved quantities (baryon 
and lepton numbers and charge), then 
three different chemical potentials are 
required. In order to get thermal and 
chemical equilibrium with the chemical 
potentials, we solve two phase balance under the Gibbs condition. 

In Fig.l, we show the boundary of the liquid-gas coexistence region of supernova 
matter and the adiabatic paths at Yl = 0.35. First, one finds that the coexistence 
region extends up to Tc ~ 15 MeV at maximum, and remains up to Tc ~ 1 MeV 
even at very low densities. These critical temperatures are much higher than those 
in asymmetric nuclear matter where Yp is fixed, and comparable to those in the 
symmetric nuclear matter. This large Tc is due to leptons: Since Yl is kept constant 
but Yp is free from constraint, the supernova matter searches its minimum free energy 
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Fig. 1. The boundary of the liquid-gas coex- 
istence region of supernova matter (thick 
solid hue), and the adiabatic paths with 
coexistence (thin solid lines) and without 
coexistence (uniform matter, dotted lines). 
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Second, all the adabatic paths seem to go across the boundary of the coexistence 
region. A Hydrodynamical model calculation shows that materials with S/B > 10 
are to be ejected to outer space. The EoS in this hydrodynamical calculation is given 
by using the TMl interaction. Therefore, the present model calculation suggests that 
high entropy supernova matter would be ejected after going through the coexistence 
phase, provided that the freeze-out temperature of supernova matter is as low as 
T = (1 — 2) MeV at densities pB = (lO"^'^ — 10~^)fm~^ and neutrinos are still 
trapped at these densities. Although the first requirement may be satisfied, the 
densities of neutrino trapping in the time scale of supernova explosion are estimated 
to he pb > lO^^fm"^. From this point of view, it would be interesting to find 
that the adiabatic path with S/B ~ 10 crosses this boundary at around {pB,T) = 
(10~^ fm~^, 1.5MeV), which satisfies both of the requirements. 

Above observations tell us that it would be marginal but still possible for super- 
nova matter experienced the liquid-gas phase transition to be ejected to the outer 
space. However, in the present model, the treatment of fragments is not realistic, 
since we have assumed that the coexisting two phases are infinite. Nuclear matter 
in the liquid phase is more symmetric than in the gas phase then both of these are 
charged, then we cannot neglect the Coulomb energy. Therefore, in the next section, 
we evaluate the fragment yield in a fragment-based statistical model. 

§3. Statistical Model of Fragments 

In order to deal with the distributions of nuclei having finite size, we choose a 
statistical model used in heavy-ion collision study, which is almost the same as the 
Nuclear Statistical Equilibrium (NSE) used in the field of astrophysics. In the calcu- 
lation, we solve the statistical equilibrium among nucleons, fragments, and leptons 
with fragment-based grand canonical ensemble. The distribution function is given 
as follows. 

Pf{A, Z) = g{T) J exp (- {Ej - pj) /T) , (3-1) 

2 

Ef = i^-Bf{A,Z) + VM.Z), pf = ApB + Zp,, (3-2) 

where Ej and Bj denote the energy and the binding energy of each nucleus, and Vc is 
the Coulomb potential energy between fragments. The spin degeneracy factor g(T) 
contains the contributions from ground state as well as from excited states estimated 
in the Bethe formula. 

In heavy-ion collision studies, the experimental values are used for the binding 
energies and only observed nuclei are considered as fragment species. In supernova 
matter, however, several modificiations are required. First, since the electron density 
is much higher in supernova matter, the Coulomb energy is screened and reduced. 
We here assume that electrons are uniformely distributed in a sphere with radius Re 
which is determined to cancel the charge of the nucleus at a given electron density. 
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Then the Coulomb energy correction is given as 

AVM, Z) = -ae^ (^^v - Iv') , V = RA/Re oc pI'^ . (3-3) 

1/3 

Because of the functional form pe , this correction is meaningfully large even at 
Pb = W~^po- For example, when the electron-baryon ratio isY^ = 0.3 the corrections 
amount to 90 % of the total Coulomb energy and 10 MeV at pB = Po and pb = 
10~^po, respectively. Therefore we cannot neglect this Coulomb energy correction 
by electrons (i.e. screening by electron) even at low densities. In addition, because 
of the finite gas nucleon density, nuclei outside of the dripline might be stabilized. 
Then we have adopted the mass table of Myers and Swiateckiu^ which contains 9000 
kinds of nuclei. 

By using the extended statistical model which uses the mass table of Ref. ^) and 
includes the Coulomb energy correction, first we have estimated the critical temper- 
atures at several densities, as shown in Table I, from the kink in the caloric curve 
and the power law fit of the fragment mass distribution. If the density is not very 
small, the freeze-out temperature is expected to be around or just below the critical 
temperature. Below the critical temperature, fragments are formed abundantly and 
nuclear number density becomes much smaller, then the mean free path for each 
nucleus becomes much longer. This rapid fragment formation makes the typical 
interaction intervals longer, and is expected to help the system to freeze-out. 
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Fig. 2. Mass distribution in supernova matter. Fig. 3. Isotope distribution in supernova matter. 

In Fig. 2, we show the fragment mass distributions at temperatures just below Tc 
for several baryon densities. When the freeze-out density is as low as pB = 10~''fm~^, 
nucleons form nuclei around the iron peak and a small number of the first peak nuclei 
of the r-process at equilibrium. As the freeze-out density increases, heavier nuclei 
are formed more abundantly. For example, nuclei around the first, second and the 
third peak of the r-process can be formed at freeze-out densities of 10~^, 10~^, and 
10~^fm~^, respectively. 

It is interesting to find that, at high freeze-out densities, the r-process peak 
structure is not clear and shifted to lower mass direction. These broadness and 
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Table I. Temperatures and typical formed nuclei shown in Fig. 2 for several baryon densities. The 
critical temperatures are written in the bracket. 



Baryon Density (fm 


Temperature (Tc) (MeV) 


Formed Nuclei 




7.0 ( 7.2) 


3rd Peak 


10"^ 


3.8 ( 4.0) 


2nd Peak 


10-^ 


1.44 (1.48) 


1st Peak 


10-^ 


0.80 (0.84) 


Fe,Ni + 1st Peak 



shifts can be well understood in the isotope distribution. In Fig. 3, we show the 
calculated even Z isotope distribution from Ca to U in the present statistical model 
at the highest freeze-out density, ps = 10~^fm~^, in comparison with the solar 
abundance. One overall normalization factor is multiplied to get a good global fit. 
We find that calculated values fit the observed (long-lived) isotope distribution very 
well from Ca to U. We also find that there are a huge number of very neutron rich 
nuclei appears in the calculated results. In the r-process, neutron absorption process 
stalls at the neutron magic numbers, but in high density supernova matter, nuclei 
beyond the drip line is statistically stabilized because of the large neutron chemical 
potential. Similar trends can be seen in the isotope distribution at lower freeze-out 
densities, although the explained range of Z is limited to Z < 56 (Ba), 40 (Zr), and 
32 (Ge) for pB = 10--\W-^, and IQ-'^fm^^. After freeze-out, these very neutron 
rich nuclei will provide a lot of neutrons by neutron emission. 

§4. Summary and Discussion 

In this paper, we have proposed the LG process which may play an important 
role for the understanding of mass and isotope distribution in the unverse. In order 
to evaluate its possibility qualitatively, we have demonstrated that the coexistence 
phase extends to low densities keeping the critical temperatures to be greater than 
1 MeV, then materials in the coexistence region can be ejected to the outer space 
at supernova within a two phase model of RMF. Next, we have performed a statis- 
tical model calculation which includes the mass table extension and the Coulomb 
energy correction to adapt to supernova matter. The calculated results show that 
nuclei above the iron peak would be formed in the coexistence region at relatively 
high freeze-out densities, and at these densities, observed isotope distribution are 
explained well. 

Works for a more quantitative discussion are in progress. 
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